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Magnetic skyrmions have recently attracted great attentions. However they are harbored in very 
limited numbers of magnets up to now. The search of new helimagnetic materials is thus an urgent 
topic in the field of skyrmion physics. In this letter, we provide a guideline on this issue, and discuss 
the possibility of realizing skyrmions in a new family of molybdenum nitrides ^ 2 Mo 3 N (H=Fe, Co, 
and Rh). By means of the first-principles calculations, the electronic and magnetic structures are 
calculated and the existence of strong Dzyaloshinskii-Moriya interaction is demonstrated. 
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The magnetic skyrmion is a swirling-like spin texture 
with nontrivial topology, where magnetic moments point 
in all directions in the space, in contrast to any other 
trivial textures such as ferrom^netism. Soon after its 
first encounter with magnetism [l| , the skyrmion was ob¬ 
served in MnSi by neutron scattering Q, and later con¬ 
firmed in Fei-xCoxSi by Lorentz transmission electron 
microscopy Q. It has attracted great attentions due to 
its promise of future applications in memory devices Q. 
Recent developments have witnessed skyrmions in several 
other bulk materials, such as FeGe hosting skyrmions 
quite close to room temperatures, and Cu 20 Se 03 Q, 
where skyrmions are insulating and exhibiting multifer- 
roic properties. 

However all these materials belong to the same lat¬ 
tice class called B20 compounds in the Strukturbericht 
Symbol. Although this class belongs to the cubic crys¬ 
tal system, complicated distributions of atoms in each 
unit cell dramatically bring down the symmetry. Thus 
the point group of B20 compounds is the tetartoidal T23 
group, the one with lowest symmetry in the cubic system, 
where both inversion and mirror symmetries are missing. 
It has been well understood that the inversion symme¬ 
try breakinggenerates the Dzyaloshinskii-Moriya (DM) 
interaction [7|, y, which competes with the Heisenberg 
exchange, and induces skyrmions or conical state at fi¬ 
nite magnetic fields. This competition is also the origin 
of the spin helices under low magnetic fields Q. That 
is why B20 compounds are termed as helimagnets. The 
DM interaction is noticeably large in B20 compounds, so 
that the dipolar interaction is negligible. This provides 
skyrmions therein with the same chirality and control¬ 
lable properties. Facing the limited choices of B20 com¬ 
pounds, an urgent problem thus arises that how to find 
other helimagnet materials harboring skyrmions, spin he¬ 
lices, or conical states driven by the DM interaction. 

In B20 helimagnets, the essential Hamiltonian in the 


continuum limit is given by 

FI = J(VS)2-f DS • (V X S) -h- S (1) 

It has been extensively tested this Hamiltonian well char¬ 
acterizes the phase diagrams of the B20 compounds. The 
last term in Eq. 0 is the Zeeman coupling. At a fi¬ 
nite window of the magnetic field and temperature, the 
skyrmion would appear, while the helix is energetically 
favored at low magnetic fields. The presence of the he¬ 
lix can be readily understood by visiting the leading two 
terms, the Heisenberg exchange and the DM interaction 
respectively, in Eq. Q- One can perform the Eourier 
transformation of the Hamiltonian, and get a quadratic 
function in momentum k from the Heisenberg exchange, 
while receiving a linear function in momentum from the 
DM interaction. The Hamiltonian is thus minimized at 
a finite momentum kg, which is the wavevector of spin 
helix. It shows a linear term in momentum is essential. 

In details, the Eourier component of the Heisenberg ex¬ 
change is given by iFHei(k;) = Jk^|Skp, where Sk is the 
Fourier component of the spin Sk = -^ / drS{r) exp(ik • 
r). Under rotations, S transforms in the same way as 
k. As a result, iFHei(k) is rotationally invariant. Fur¬ 
thermore, inversion and mirror symmetries are also re¬ 
spected. Therefore this is a generic quadratic term for 
all ferromagnets. 

On the other hand, the DM interaction provides a lin¬ 
ear term in momentum; iFDM(k;) = iDSk ■ (k x S_k)- 
Although the rotational symmetry is still preserved, in¬ 
version symmetry is apparently broken, which is the well 
know precondition for the DM interactions. However a 
long overlooked fact is the mirror symmetry is also bro¬ 
ken by this DM interaction. This comes from the fact 
that Sk • (k X S_k) is a pseudoscalar. Under any im¬ 
proper rotation such as mirror reflection in the lattice, 
this term flips sign, and should be ruled out in the en¬ 
ergy. That is why magic is witnessed in B20 compounds. 
The point group T23 only has pure rotations, thus al- 
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lows the DM interaction in this form. In case a mirror 
plane exists, the DM interaction arise from broken inver¬ 
sion symmetry must be staggered, failing to end up with 
a continuum limit in this form. Now a question arises 
whether this form is the only allowed term linear in k for 
any materials. 

To answer this question, we can rewrite any k-linear 
terms as a tensor product iJDM(k) = idijmk^ , 
where dijm is a third order tensor that can be constructed 
from symmetry analysis. Any symmetry operation R can 
be represented as a 3x3 matrix in natural basis (x, y, z). 
Under such operation, vector k transforms as ki —> kjRji, 
while the pseudovectors S±k transform as Si —>■ \R\SjRji, 
where |i?| is the determinant of R matrix. If R is an im¬ 
proper rotation, |i?| = —1. Once i? is a symmetry oper¬ 
ation, energy should be invariant under such rotation, 
therefore the tensor dijm must satisfy the Neumann’s 
principle: 

dijm — ddipRjqRmrdpqr ( 2 ) 

In practice, one does not need to go through all symme¬ 
try operations in order to determine the d tensor. Most 
operations can be written as products of some indepen¬ 
dent matrices, called generating matrices 0 , within the 
same point group. For T23 point group, the generating 
matrices are C 2 and C 3 rotations. The Neumann’s prin¬ 
ciple thus leads to the constrain that d^yz = dyzx = d^xy, 
and dxzy = dyxz = dzyx- One can symmetrize these pa¬ 
rameters by dxyz = S + D and dxzy = S — D. How¬ 
ever because the whole Hamiltonian can be reorganized 
as I]k^DM(k) = Y^kKdijm - dimj)k^Si^Sl^, the sym¬ 
metric component S does not contribute. The resulting 
Hamiltonian is thus '^^^iDeijmk'^ , which repro¬ 
duces the DM interaction in Eq. Q- 

The same method applies to any other lattices. Con¬ 
tribution to the Hamiltonian from any tensor with re¬ 
dundant indices vanishes when completing the summa¬ 
tion over momenta. The relevant terms are six compo¬ 
nents with indices permutations of (x, y, z). For future 
convenience, these six components are symmetrized as 

dxyz — “b dyxz — dyzx — “b 

dxzy = Ps - /3a, dzxy = 7S + 7A, dzyx =75“ 7A• As a 
result, the total Hamiltonian is given by 


u 

CO 

II 

[J(VS)2-h-S 

-b 

i(as -b OA - /3s + /3 a)S • {dx x S) 

+ 

^(-as +aA + ds + I3 a)S ■ {dy x S) 

-b 

7 aS • {dz X S)] 


where dr is the directional derivative along r direction. 
Under low magnetic field h, spin helix is thus formed 
along certain directions given the competition between 
anisotropic DM interaction and the Heisenberg exchange. 


TABLE I: Constrains of nonzero dijk parameters for all pos¬ 
sible point groups 


Class 

Constrains 

Point Groups 

I 

No Constrain 

Cl, C 2 , D 2 

II 

as = /3s = 7 s = 0 

C 4 , D 4 , C 3 , D 3 , Ce, De 

III 

aA = Pa = yA = 0 

Si, D 2 d 

IV 

ois = Ps = ys, ctA = Pa = 7a 

T 

V 

as = Ps = ys = 0, OA = Pa = 7 a 

0 


A complete list of point groups contributing to nonzero 
DM interactions is summarized in Table HI where the B20 
compounds are located in class-IV. We know that the 
point group can be decomposed into two categories, the 
ones with and without improper rotations. Most groups 
in Table [J belong to the category without any improper 
rotations. Mirror, in addition to inversion, is broken in 
these lattices. The only exception is class-HI; S 4 and 
D 2 d groups, where S • (9z x S) is prohibited. Helices 
or skyrmions can only form in the plane perpendicular 
within the horizontal mirror. In contrast, class-H allows 
the presence of the DM interaction in the z-direction. 
The DM interaction of this class can be rewritten as 
i(aA+/ 3 A+ 7 A)S-(VxS)-f i 7 AS-(i 92 xS). It shows explic¬ 
itly that the z-direction is distinct from the xy plane. Al¬ 
though helices or conics can propagate along z-direction, 
the skyrmion crystal would favor to accommodate in the 
xy plane instead. However one should be aware that the 
spin anisotropies are quite large in these lattices with 
reduced symmetry, which are enemies of the skyrmion. 
One can also construct spin anisotropies by the symme¬ 
try tensors up to arbitrary order. These details are out of 
the scope of this work. In the future, case by case studies 
can be conducted for specific materials. 

The most important message delivered from Table U is 
class-V, which is described by exactly the same Hamil¬ 
tonian as the B20 compounds in class-IV. The spin 
anisotropies are also the same in these two classes. There¬ 
fore we expect the spin physics observed in B 20 com¬ 
pounds are also persistent in the point group O. There 
exists 8 space groups in O group. We found that the 
most promising material in analogy to B20 compounds 
is the A 2 M 03 N family with A=Fe, Co, Rh, or their alloys. 
The pure A 2 M 03 N 13-1^, shown in Fig. [HJa), has the 
filled /3-manganese structure, where A atoms lie on the 
8 c positions of a cubic unit cell with space group symme- 
try P4i32 forming a single (10, 3)-a network jl6| . The 
space within this network is filled by corner-shared MoeN 
octahedra. This structure is entirely analogous to that 
of /3-manganese, differing only by the addition of the in¬ 
terstitial non-metal atoms. We expect large strength of 
spin-orbital coupling in this family, and thus skyrmions 
can be hosted. 

To quantitatively confirm the aforementioned conjec¬ 
ture of promising material, we carried out the first- 
principle calculations of pure A 2 M 03 N using the pro- 
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FIG. 1: (Color online) (a) The schematic crystal structure of 
the filled /3-manganese structural A 2 M 03 N (T=Co, Fe, and 
Rh). Space within the (10, 3)-a network of A atoms (green) 
is filled by vertex-sharing MoeN. (b) The charge density dis¬ 
tribution in the (001) plane crossing the A-}Ao-A atoms in the 
NM state calculations. 



Energy (eV) 


jected augmented wave method as implemented in the 
VASP code [m , where the exchange-correlation potential 
was calculated using the generalized gradient approxima¬ 
tion (GGA) as proposed by Pedrew, Burke, and Ernzer- 
hof (PBE) [l 2 |. All atomic positions and lattice constants 
of A2M03N were allowed to relax simultaneously to mini¬ 
mize the energy. A 500 eV cutoff in the plane wave expan¬ 
sion ensures the convergence of calculations up to 10“^ 
eV, and all atomic positions and the lattice constants 
were optimized until the largest force on any atom was 
0.005 eV /A. Eurthermore, we used a 8 x 8 x 8 Monkhorst- 
Pack k-grid Brillouin zone sampling throughout all of cal¬ 
culations. In addition, the spin-orbit coupling was also 
included with the second variational method. The lat¬ 
tice constants and the internal coordinates of the atomic 
positions for the systems of A2M03N are all optimized 
and listed in Table IE which shows that the optimized 
lattice parameters for both Ee2Mo3N and C 02M03N are 
quite consistent with that from experiments |13l4l5l |. Al¬ 
though lattice parameters for Rh2Mo3N are lacking ex¬ 
perimentally, our optimized values are larger than those 
of Pe2Mo3N and G02M03N, which are reasonable as the 
atomic radii of Rh is apparently larger than that for Ee 
and Go. 

The electronic properties of A2M03N in the quenched 
paramagnetic state in A 3 ( 4 )d and Mo Ad orbitals are 
studied. Such studies provide references for forthcoming 
magnetism studies. By analyzing the density of states 


TABLE II: The optimized lattice constants as well as the 
internal coordinates for the filled /3-manganese structured 
A 2 M 03 N (A=Co, Fe, and Rh). 


^2Mo3N 

A=Fe 

0 

0 

II 

A=Rh 

a (A) 

6.6458 

6.6367 

6.8224 

A{x = y = z) 

0.0723 

0.0677 

0.0598 

Mo {y) 

0.2019 

0.2002 

0.2021 

Mo {z) 

0.4519 

0.4502 

0.4521 


FIG. 2: (color online) Total DOS and PDOS on A 3(4)d, Mo 
Ad, and N 2p orbitals of the NM state for A 2 M 03 N with (a) 
A=Fe, (b) A=Co, and (c) A=Rh. The Fermi energies are set 
to zero. 


(DOS) at the Fermi level we can infer whether the mag¬ 
netic state is favored. Fig. O shows the total DOS of 
A2M03N, and its projected DOS (PDOS) onto A 3 ( 4 )d, 
Mo Ad, and N 2p orbitals respectively. Furthermore, it 
also shows that the mixing between A 3 ( 4 )d and Mo Ad 
occurs mainly around the Fermi energy ranging from —4 
eV to 2 eV indicating the sizable 3(A)d-Ad hybridization 
between A 3(A)d and Mo Ad orbitals. Particularly, com¬ 
parable contributions to the conducting carriers are wit¬ 
nessed from the Rh Ad and Mo Ad orbitals in Rh2Mo3N, 
shown in Fig. [ 21 (c). This is attributed to the much more 
extended Ad orbitals of Rh compared to the 3d orbitals 
of Fe and Go, which leads to a stronger hybridization, 
expanded bandwith, and highly overlapping between Rh 
Ad and Mo Ad orbital states. Additionally, it is impor¬ 
tant to point out that the DOS coming from A 3(A)d 
orbitals is nonvanishing at the Fermi level. The val¬ 
ues of DOS at the Fermi level are ~ 

= 0 . 73 ], lNS:(Ef) = 2.3 and = 0 . 8 ], 

and = 1.5 and N^^{Ef) = 0 . 67 ] states 

per eV per A(Mo) atom for Fe2Mo3N, G02M03N, and 
Rh2Mo.3N, respectively. According to the Stoner crite¬ 
rion [ITI, [l8j, magnetism may occur only when the DOS 
satisfies N{Ef)I > 1 , where I is the Stoner parameter, 
which takes values of 0.7 eV- 0.9 eV for ions near the mid¬ 
dle of the transition metal series (note that the effective 
/ can be reduced by hybridization). It shows clearly that 
although the NM state is favorable on the Mo atoms, it 
is unstable against the magnetic states on the A atoms. 
A atoms are responsible for magnetisms in this family. 

To explore the magnetic behaviors of A atoms in 
A2M03N, we consider a local exchange model based on 
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the nearest neighbor Heiseirberg and DM interactions: 

H = J2 • S, + E • (S. X S,), (4) 

{i.J) (i.J) 


where is the operator of A spin at site i, (i,j) de¬ 
notes the summation over the nearest neighboring sites 
between A atoms. Parameters Jij arrd are the nearest 
neighbor Heisenberg and DM interactions, respectively, 
which can be evaluated by using the four-state energy- 
luapping analysis 0. The results are listed in Table Hill 
Here lists only the coupling constants oir the three bonds 
shown in Fig. [TJa). Interactions on other bonds are the 
same except for a rotation. The vector of the DM 
interaction on each bond are almost perpendicular to the 
bond as expected. That is because the lack of inversion 
symmetry generates a local electric field e on each bond. 
An electron hopping between two ends of this bond me¬ 
diates the exchange interaction of the neighboring spin, 
and feels an effective magnetic field b = v x e, where 
velocity v is along the bond. Spin of this electron thus 
proceed about b, resulting in a DM interaction of neigh¬ 
boring spins with the DM vector parallel b. Therefore 
Dy should be perpendicular to the bond. The small 
deviation from perpendicular is associated with the non¬ 
uniformity of e. The consistency between this physical 
picture and the calculation results justifies our evalua¬ 
tion of the DM interactions. From these data we no¬ 
tice that the Heisenberg interactions Jij are almost the 
same on the three bonds in each compound. They reach 
the maximal value in Rh 2 Mo 3 N, which attributes to the 
nature of strong hybridization between Rh Ad and Mo 
Ad orbitals [see the charge density distribution shown in 
Fig. mb)]. The ratio between the DM interaction and 
Heisenberg exchange has the largest value 7 = |-j^| ~ 
0.084 in C 02 M 03 N, while that of Fe 2 Mo 3 N and Rh 2 Mo 3 N 
are 0.069 and 0.073, respectively. These values are much 
larger than conventional ratios of 7 < 0.05 Such 
strong DM interaction between A-A atoms in A 2 M 03 N 
mainly comes from strong spin-orbital coupling in Mo Ad 
orbitals, which mediate the spin-spin interaction between 
A atoms as indicated from the charge density distribution 
shown in Fig.[T][b)- Thus, we conclude that the C 02 M 03 N 
with strong DM interaction is a promising candidate for 
realizing the exotic skyrmion in pure A 2 M 03 N series. 

As the skyrmion radius is controlled by the ratio 
between the DM interaction and Heisenberg exchange, 
it is of great interesting, from both physics and ap¬ 
plications’ perspectives, to find materials with large 
7 = To this end, material optimizing is re¬ 

quired. Here, we suggest that the substituted compound 
Co 2 -a;-j/Rha;FeyMo 3 N has a larger 7 value than pure 
C 02 M 03 N by fine tuning the stoichiometry. According 
to the aforementioned calculations and discussions, the 
strength of DM can be raised through an enhancement 
of the hybridization between A 3(4)d and Mo Ad orbitals. 


TABLE III: The values of magnetic exchange coupling and 
DM interactions on the three bonds shown in Fig. [Ha) in the 
filled d-manganese structured A2M03N (A=Co, Fe, and Rh). 
The unit is meV/S^, where S is the spin of the A atom. 



II 


0 

0 

11 

2=1 

Di 

II 

Jl 

-27.17 


-41.46 

-62.71 

J 2 

-27.23 


-41.46 

-62.85 

J 3 

-27.25 


-41.46 

-62.84 

Di 

( 0.36,-0.37, 

1.80) 

( 1.19,-1.16, 3.06) 

( 2.39,-2.50, 1.71) 

D 2 

^0.53, 0.61,- 

0.07) 

(-1.05, 1.12, 2.69) 

^1.51, 1.26, 4.35) 

Da 

( 0.61, 0.48, 

0.08) 

( 1.11, 1.05,-2.69) 

( 2.02, 1.41,-4.31) 


Therefore partially substituting Co by Rh will help as the 
Ad orbital of Rh is more extended than the Zd orbital of 
Co, which makes the overlapping to Mo Ad orbital more 
sufficient. However, such substitution simultaneously in¬ 
troduces large atomic orbital potential difference between 
large radii Rh and small radii Co atoms, which gener¬ 
ates large scattering on itinerant electrons and weakens 
the exchange interactions. To overcome those potential 
barriers, an isovalent dopant Fe needs to be introduced. 
Based on this analysis we further calculate the magnetic 
exchange interactions for CoRho, 75 Feo. 25 Mo 3 N and ob¬ 
tain the ratio 7 = up to 0.11. Therefore, we suggest 
the substituted compound Co 2 _a,_yRh 2 ,FeyMo 3 N is the 
most promising material realizing the exotic skyrmion 
state. An experimental hint has already been observed 
in [isl , where zero remenance and a kink of magnetic sus¬ 
ceptibility both indicate the presence of spin helices, the 
forerunner of skyrmions, at small magnetic fields and low 
temperatures. A complete exploration of the phase dia¬ 
gram and stoichiometry is required. On the other hand, 
the Curie temperature can be elevated by introducing 
other dopants such as Pt 14|. 


In conclusion, we have constructed a framework of 
searching new helimagnet materials harboring exotic spin 
textures of skyrmions. The effective Hamiltonian is 
derived based on symmetry analysis. A new family 
A 2 M 03 N (A=Fe, Co, and Rh) has been proposed, and 
the first-principle calculations are performed. We hope 
that this family is just a corner of a huge iceberg. 
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